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ABSTRACT: This study describes a generalized analog computer program 
for control of the feed stream enthalpy to a distillation column. It is based 
on the solution of a steady state heat balance equation written for the feed 
system to the column. Control is accomplished by feeding the solution 
from the computer to a conventional recording-controlling instrument 
which automatically manipulates a control valve in the stream line to a 
feed preheater. 

The program developed is designed for solution with an EAI PC- 12 
Process Control Computer. The programming flexibility of this computer 
permits the use of standard, solid-state computing components. Scaling 
of the circuit is completely general, enabling the program to be utilized 
for a variety of columns having different operating conditions. 



GENERAL 

The importance of good regulation and control of 
the three major heat inputs to distillation columns 
in achieving improved operations is cited in the 
literature (1,2). Among the reasons presented 
therein are: 

(1) the difficulty of separation in many columns 

(2) the change in the dynamic character of 
columns 

(3) the interaction of variables internal to the 
column, or external due to column auxili- 
aries 

(4) the non-linear nature of significant varia- 
bles in a column 

(5) the change in demand upon column oper- 
ation due to its dependency upon other 
parts of the process. 

Even though control of distillation columns is a 
well established art, and is adequately docu- 
mented (3,4), there are certain heat input dis- 
turbances which require more than conventional 
instrumentation for their control. One such dis- 
turbance is the change in enthalpy— heat content— 
of feed entering the column. 

COMPUTER CONTROL OF FEED ENTHALPY 

Although the heat supplied to a column by the 
feed is usually small in comparison to the total 
heat required for the separation, the importance 
to good column control of its being stabilized 
should not be overlooked. It is generally desirable 



that the feed enthalpy be controlled at a level 
which will result in the least operating costs 
while producing on- specification products. Re- 
gardless of whether the feed must enter at its 
bubble point temperature, partially vaporized, 
or subcooled, the problem of regulating it at 
a selected value becomes important to good 
column control. 

One practical and economical method for regu- 
lating the feed heat content involves a special 
purpose analog computer, called a Feed Enthalpy 
Computer. This approach has been found, by 
simulation studies and extensive plant tests (1), 
to give excellent results where the feed enters 
at its bubble point temperature, partially 
vaporized. 




Figure 1. Feed Enthalpy Computer Control 
of Distillation Column. 
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The concept of the Feed Enthalpy Computer is 
based on the inferential measurement of the total 
feed heat content by adding appropriate heat 
quantities; i.e., the initial heat content of the 
feed entering the economizer exchanger, the heat 
supplied to the feed by the economizer exchanger, 
and, the heat supplied by the steam feed pre- 
heater (refer to Figure 1). The computer con- 
tinuously calculates the feed heat content in BTU 
per pound of feed at the column input. Its 
output serves as the measured variable to an 
enthalpy recorder-controller whose output cas- 
cades into a steam flow controller that modulates 
the preheater steam flow to maintain a constant 
feed enthalpy. Inputs to the computer are pro- 
vided by conventional instrumentation in the form 
of signals from flow transmitters and thermo- 
couples. 

Benefits which accrue from improved column 
control through regulation of the feed enthalpy 
by computer control are (1): 

(1) Improved primary control of the separation 

(2) Lower operating costs 

(3) Higher column throughputs 

(4) Smoother terminal stream flows 



balance equation written for the column feed 
system. Referring to the schematic diagram of 
the feed portion of the distillation column (Figure 
1) the total feed heat content per pound of feed, 
neglecting losses, is the sum of the following 
heat quantities. 
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= Initial feed heat content above 
some reference temperature T Q 
— BTU per pound of feed. 

T B ) = The heat given up by the 
2 bottom product stream to 
the feed in the economizer 
exchanger — BTU per pound 
of feed. 



H 3 = — H g = The heat given u p to the feed by the 
1 steam at the feed preheater — BTU per 

pound of feed. 

The total heat content of the feed stream, per pound of 
feed, is then, in equation form; 



H T= H 1+ H 2+ H 3 



(1) 



(T 1 -T o ) + 
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<T R -T R ) + — H a (2) 



ANALOG COMPUTER CIRCUIT 



(5) A more predictable operation with less 
human attention. 

DERIVING THE EQUATION 

Control of the feed stream heat content is based 
upon the computer solution of a steady state heat 



A suggested solution procedure outlining the 
mathematical operations required for solving 
equation 2 is shown in the information flow sheet 
of Figure 2. Using this diagram it is a relatively 
simple task to construct the computer circuit 
of Figure 3. Note that this circuit assumes that 
signals in the form of dc currents are available 
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Figure 2. Calculation Scheme for Solving Steady-State Feed Heat Content Equation with an 

Analog Computer. 
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Figure 3. Scaled Analog Computer Circuit for Solving Feed Heat Content Equation, i. e. , General Purpose Program for 

Feed Enthalpy Computer. 



for representing the square of the three flow 
rates. Square root operations on these flow signals 
are performed automatically, being programmed 
into the computer circuits. Temperatures are 
assumed to be measured by thermocouples whose 
outputs may be fed directly to amplifiers internal 
to the computer. Although the computer output as 
shown in Figure 3 takes the form of a dc voltage 
which varies from 0 to -10 volts it is possible 
to arrange the amplifier input and output circuits 
such that a dc current acceptable by conventional 
instruments is obtained. 

The circuit of Figure 3 employs the following types 
of standard PC- 12 Computing Components. 

DUAL OPERATIONAL AMPLIFIER (Type 6.368). 
Provides two high gain, chopper stabilized, tran- 
sistorized operational amplifiers for the per- 
formance of linear and non- linear mathematical 
operations. Low drift, high accuracy, field-proven 
circuits are designed specifically for general 
purpose on-line computing applications. Compact, 
plug-in modules mount in vapor-tight industrial 
housings. Used in the computer circuit of Figure 3 
for scaling of input signals and for summing 
operations. Four dual modules or a total of eight 
amplifiers are required. 

LOW LEVEL DIFFERENTIAL AMPLIFIER (Type 
6.422). A transistorized, low-level differential 
amplifier designed for amplification of low-level 
signals. Design features include AC input section 
for added stability, 0.01% gain accuracy, and input 
impedances greater than 10 megohms. Used in the 
computer circuit of Figure 3 to accept and amplify 
differential thermocouple inputs for measurement 
of the bottoms flow differential temperature. One 
required. 

THERMOCOUPLE AMPLIFIER (Type 6.406). 
Transistorized, low-level amplifier designed for 
amplification of low level thermocouple signals. 
Has similar characteristics as Type 6.422 Low 
Level Differential Amplifier but is designed speci- 
fically for single-ended input signals. Used in the 
circuit to accept and amplify a thermocouple 
input for measurement of the input feed temper- 
ature. One required. 

COEFFICIENT SETTING POTENTIOMETERS 
(Type 12. 779). Provides four 10-turn, wire wound, 
2000 ohm potentiometers with calibrated dials for 
setting of equation constants, constant inputs, and 
bias voltages. Six potentiometers required. Used 
in the circuit for introduction of scaling constants 
and for setting values for T Q , H g , Cp^, Cp^. 

1/2 LOG X DIODE FUNCTION GENERATOR 
(Type 16.215). Solid-state dual fixed diode func- 
tion generator for generating logarithms, ex- 
ponential functions, etc. When used with Ampli- 
fier 6.368 one generator will accept a positive 
voltage and generate a straight- line approximation 
of the logarithmic curve. The other generator 



accepts a negative voltage to perform the same 
function. Internal circuit components are sized 
to produce an output of 2.5 Logio(lOX) for an 
input of X to facilitate the combined operations 
of square root and multiplication with a single 
amplifier. Current output for an input of X is 
5 x 10~ 4 Log-iQ(10X). Used in the circuit for form- 
ing the log of the signals representing the square 
of the flows. Three negative input generators 
are required. 

LOG X DIODE FUNCTION GENERATOR (Type 
16.214). Similar to the 1/2 Log X Generator but 
with internal circuit components sized to produce 
an output of 5.0 Logio(lOX) when used with an 
operational amplifier. Current output for an input 
of X is 10" 3 Logio(lOX). Used in the circuit to 
perform anti-log operations and to take the log 
of the signal representing the bottoms differential 
temperature. Two positive input and one negative 
input generators are required. 

These components, together with appropriate am- 
plifier input and feedback resistances, system 
power supplies and control modules, patching sys- 
tem, etc., are mounted in an industrial housing 
to form the Feed Enthalpy Computer. 

COMPUTER CIRCUIT SCALING 

A generalized scaling scheme has been devel- 
oped for the computer circuit used to solve the 
steady- state heat equation. All signals and circuit 
components are expressed in terms of the ex- 
pected maximum values of the input variables. 
Thus, once the column operating conditions are 
known, the circuit components and the computed 
signals can be completely specified using the 
scaled computer variables shown on the circuit 
of Figure 3 and a design equation developed from 
the scaling procedure. 

The Appendix presents a detailed development of 
the scaling. 

EXAMPLE 

Assume that the column operating conditions are 
given as follows: 



AT T 



= 250,000 pounds per hour 
= 300°F 

= 500,000 pounds per hour 

= 20,000 pounds per hour 
= 200°F 

= 0.532 BTU per pound per °F 
= 0.424 BTU per pound per °F 
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= 750 BTU per pound of steam 



#5 _ T q /T 1 75/200 - 0,375 



T Q =75°F 

From equation 2, the maximum value of is calculated as 
250,000 

H rj, = 0.532 (200) + (0. 424X300) + 



20,000 
500,000 



500,000 



(750) = 106.4 + 63.6 + 30 = 200.0 



#6 - c = 0,532 
P F 



Shunt resistors R g are calculated to have the value . 
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Through use of the design equation developed in the Appen- 
dix .... 



F ma x K 2 



F max K3 



10 



K c AT R l n 3 s 

5 B max B max iU ^max 



Then 
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In addition 



AT 



20 F„ 



max B max = (300X250,000) _ 7 5 



(20X500,000) 
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= 625 ohms 



Biasing resistors R are calculated to have the value .. 



R = 



10 



R s l B 



-(10,000) = 



10 4 



(625)(4 x 10 - 6 ) 



40,000 ohms 



Scaling resistor has the value 



10 K c 
1 = ~, 



^5 (10)(20) 

R,= (10,000) = —— (10,000) 

200 
A max 

= 10,000 ohms 

Amplifiers of interest have the following scaled com- 
puter variables as outputs 



and 



H 



K, 



10 3 S 



max = (1000X20,000) 2 Q 



20 F max (20X500,000) 



Potentiometer Settings become 

#1 - 0.5 Log 10 (10 c ) = 0.5 Log 1Q (4.24) = 

(0.5X.6274) = 0.3137 
#2 - 0.5 Log 10 (100/K 2 ) = 0.5 Log 10 (13.33) = 

(0.5X1.1248) 
= 0.5624 

#3 - 0.5 Log 10 (10 K 3 ) = 0.5 Log 1Q (20) = 
(0.5X1.301) = 0.625 

#4 - H / 1000 = 750/1000 = 0.750 
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. . . while potentiometers #4 and #6 have outputs 
scaled as follows 
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-NOMENCLATURE - 

temperature of feed before entering econ- 
omizer; °F. 

arbitrary reference temperature used to 
compute H T ; °F. 

temperature of bottoms product entering 
economizer; °F. 

temperature of bottoms product leaving 
economizer; °F. 

feed flow rate; lb/hour. 

bottom product flow rate; lb/hour. 

steam flow rate; lb/hour. 

difference in enthalpy of steam entering 
preheater and condensate (assumed con- 
stant); BTU/lb of steam. 

average specific heat of feed; BTU/lb — 
°F. 



average specific heat of bottom product; 
BTU/lb - °F. 



K 2 , K 3 , K 5 = scaling constants 
R 



R 



shunt resistors used to develop voltage 
from current inputs; ohms. 

biasing resistors on inputof flow scaling 
amplifier; ohms. 

input resistor for Hj input to amplifier 
#10; ohms. 

current representing square of input flows 
millamperes. 



SUBSCRIPTS 

max = expected maximum value or range of variable. 
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- APPENDIX - 



INPUT SCALING AMPLIFIERS „ The components of the following circuit must be sized such that, for the 
range of currents between the minimum and maximum, 



e = 0 for Id . 
o J^mm 
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Summing currents at amplifier grid: 



Solving for resistance R: 
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Therefore, when e Q = 0 and 1^ = I B 



mm 
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The current summation equation then becomes 



Vb 10 _fo_ _ 0 

10K 100K s B min + 10K 



so that when I B = I Bmax > % = -10V 
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e Q is to be made to represent the square of a flow, say B , so computer variable is defined as 

-e = -k B 2 
o 

where k is an amplitude scale factor. But 

k • B 2 = 10 
max 

so 



k = 



10 



B 



max 



The scaled computer variable thus becomes 



max 



The Thermocouple Amplifier and the Differential Amplifier are scaled in exactly the same manner so as 
to restrict their maximum output to 10 volts. The circuits are more complex, however. Refer to the 
appropriate Service Manual for the exact procedure. 

SCALING OF AMPLIFIER OUTPUTS 

Amplifier #6: Summing Currents at grid (See Figure 3) 
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Amplifier #7: Summing currents 
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Amplifier #8: Summing currents 
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Amplifier #9: Summing currents 
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Amplifier #10: Summing currents 
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In order to sum these signals, all scale factors must be equal, or 
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The equation then becomes 
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i^[ H l + H 2 + H 3] = - e o = ^ 
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Representing as a computer variable, and restricting the amplifier output to a 10 volts maximum 
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Thus, the design equation becomes 
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